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Abstract: 

The gluon polarization AG in a nucleon can be defined in a gauge invariant way as the 
integral over the Ioffe-time distribution of polarized gluons. We argue that for sufficiently 
regular polarized gluon distributions AG is dominated by contributions from small and 
moderate values of the Ioffe-time z < 10. As a consequence AG can be estimated with 
20% accuracy from the first two even moments of the polarized gluon distribution, and its 
behavior at small values of Bjorken x or, equivalently, at large Ioffe-times z. We employ 
this idea and compute the first two moments of the polarized gluon distribution within the 
framework of QCD sum rules. Combined with the color coherence hypothesis we obtain 
an upper limit for AG ~ 2 ± 0.5 at a typical scale /i 2 ~ 1 GeV 2 . 
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During the last two decades deep-inelastic lepton-nucleon scattering (DIS) has proven 
to be one of the most valuable sources of information about nucleon structure. According 
to the Operator-Product-Expansion (OPE), the Q 2 -dependence of the cross-section can 
be entirely included in perturbatively calculable Wilson coefficients, and thus completely 
factorized from the effects of long-distance dynamics, described by a set of twist-2 parton 
distribution functions which are defined at a certain reference scale. The analysis of 
experiments with unpolarized lepton beams has provided detailed informations on quark 
momentum distribution functions and has lead to the conclusion that - even at relatively 
low scales of few GeV 2 - the nucleon momentum is shared almost equally between quarks 
and gluons. With the advent of high-quality polarized targets and beams it has become 
possible to study polarized quark distributions 0. In particular, combining experimental 
results with SU(3)-flavor symmetry allows to determine the fraction S of the nucleon 
polarization carried by quarks @. Since the latter turns out to be much smaller, S ~ 0.3, 
than expected from simple quark models || , it is natural to search for other contributions 
to the nucleon spin. 

At a given scale fi 2 the nucleon spin sum rule allows to separate the nucleon spin into 
gauge-invariant quark and gluon contributions [|J: 

~ = J,(/i 2 ) + J,(/i 2 ) , (1) 

with both parts calculable as matrix elements of local operators sandwiched between 
polarized nucleon states. The quark contribution can be split further into gauge-invariant 
polarization and angular momentum parts: 

J q {f) = \n» 2 )+L q {n 2 ). (2) 

Here the quark polarization S(/^ 2 ) is linked to the matrix element of the flavour singlet 
axial current operator / g? 3 £'0775'0, and the angular momentum L q (^ 2 ) is connected to 
the gauge-invariant angular momentum operator — / d 3 xip(x x iD)ip, with the covariant 
derivative _D M = (Do, D) = <9 M + igN 1 . The separation of the gluon contribution J g into 
angular momentum and polarization parts, 

J a Qj?)=AG(vL 2 ) + L g (jj?), (3) 

is more involved. The reason is that in QCD appropriate relations between (even) mo- 
ments of the polarized gluon density and matrix elements of twist-2 operators exist only 
for I > 2: 

f 1 dxx 1 AG(j,/i 2 ) = TiU 2 ), with 1 = 2,4,..., 
Jo 

(p,s\ nil n u TrG^{in- D) l - X & v \p, s) = (s ■ n){p ■ n)^(/i 2 ). (4) 

Here x stands for the Bjorken scaling variable which is equal to the gluon light-cone 
momentum fraction and G^ u is the gluon field strength tensor. The nucleon state \p, s) 
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carries an invariant mass m^v, momentum p and spin s, with s 2 = —m? N . The light-like 
vector n M , n 2 = 0, projects out the "plus" component of any four vector, i.e. n ■ a = 
a + = a + a 3 . There is no corresponding relation for I = or AG(/i 2 ) = Jq 1 c&e AG(x, /i 2 ), 
because a suitable gauge invariant, charge conjugation even, local operator does not exist. 
Thus contrary to the quark case, the gluon polarization cannot be extracted from a matrix 
element of a local QCD operator. 

On the other hand experiments, like large Pt jets or open charm production in polar- 
ized DIS, can in principle provide experimental informations about the polarized gluon 
distribution AG(x, /i 2 ) (see e.g. ||) and its integrated strength AG which, as a measur- 
able quantity, should have a gauge-invariant definition. Indeed, consider the non-local 
twist-2 operator |J 

0(A/2; -A/2) = n^Tr G*(A/2)[A/2; -A/2]G£(-A/2), (5) 

where A is a light-like vector proportional to n. To ensure gauge invariance the operator 
fl5|) contains the path-ordered exponential in the adjoint representation [A/2; —A/2]. The 
nucleon matrix element of the operator (||), renormalized at a scale /i 2 , can be written in 
terms of the polarized gluon density AG(x,/i 2 ): 

(p, s\0(A, 0)\p, s) — (p • n)(s ■ n) / dx x AG(x, fi 2 ) sin [x(p ■ A)] , (6) 

Jo 

which defines the corresponding distribution T(z, fi 2 ) in terms of the "Ioffe-time" z = p- A 
0as: 

T(z,[i 2 )= / dxx AGix, ii 2 ) sinixz). (7) 
Jo 

Performing a Taylor expansion of eq. @ around A = one can easily see that the definition 
of the Ioffe-time distribution is equivalent to the OPE relation between matrix elements 
of local operators (§) and moments of AG(x, /j, 2 ). Now, although the first moment of 
AG(x, fi 2 ) is not determined by a local operator, it can be extracted under the assumption 
that AG(x,fi 2 ) is sufficiently regular at x — 0. Using the relation: 

roo i 

/ dz sin(xz) = P— , (8) 
Jo X 

we obtain: 

r\ /*oo 

AG(fi 2 ) = / dxAG(x,n 2 )= dzT{z,fi 2 ). (9) 
Jo JO 

In the the Fock-Schwinger gauge the integral ([|) can be related to the matrix element 
of the topological current K* = ffi e^P a A a u (d p A a a + \gf abc A b p A c a ) §. Recently Jaffe @ used 
this relation to compute AG(/i|) in the bag model, taking into account gluonic contribu- 
tions arising from one-gluon exchange between quarks. In a semi-classical framework, i.e. 
neglecting renormalization, one obtains a negative value of AG ~ —0.2 at the bag scale 
\j? B . In the present paper we propose another procedure which allows, under favorable cir- 
cumstances, to estimate AG(/i 2 ) with an accuracy sufficient from any phenomenological 
point of view. 
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Equation (^) defines the gluon polarization in terms of all higher moments of the 
polarized gluon distribution, and as such has a limited practical valueQ. We argue, how- 
ever, that a reasonable estimate for AG(/i 2 ) can be obtained from the first few moments 
only, assuming a sufficiently smooth polarized gluon distribution. In particular our pre- 
vious studies []nj of unpolarized Ioffe-time distributions have shown that an overall good 
approximation can be obtained by combining small and large z expansions and the be- 
havior of parton distributions at large values of the Bjorken variable x — > 1. At small 
longitudinal distances, i.e. small z, the Ioffe-time distribution can be approximated by 
its Taylor expansion around the origin, with coefficients calculable as matrix elements of 
local twist- 2 operators of the lowest dimension. At large z its shape is sensitive to the 
physics of large longitudinal distances, which corresponds to small x. In this region a 



proper description of parton distributions in terms of QCD is still under debate [11]. An 



analysis of unpolarized distributions |L0j shows that the transition between the small and 
large z domain is relatively sharp, and that the large longitudinal distance regime prevails 
beyond z ~ 10. Note that within a classical picture of a nucleon localized at the origin a 
Ioffe-time of z ~ 10 corresponds to a longitudinal distance of approximately 2 fm which, 
perhaps not surprisingly, coincides with the classical nucleon sizeQ. 

In Figure 1 and 2 we show Ioffe-time distributions resulting from two very different 
polarized gluon densities. While the distribution of Brodsky et al. [0 gives AG = 0.54 
at a typical hadronic scale, the distribution of Chiappetta et al. |H| leads to AG = 1.7 
at fj, 2 = lGeV 2 . In both cases the correlation function T(z,^ 2 ) falls off quickly at large 
z and, due to the resulting triangular shape of T(z, /x 2 ), the main contribution to the 
first moment (^) is contained in the region z < z , with zq around 10. This agrees with 
the intuition that degrees of freedom which determine the large z behavior of parton 
distributions do not couple to the target polarization. Note that if AG(x,/j, 2 ) behaves 
like ~ l/x a one finds T(z,fi 2 ) ~ l/z 2 ~ a . The dominance of contributions from moderate 
z can be quantified by calculating the fraction of AG(/i 2 ) which results from the interval 
between and a given z, as shown in Figure 1 and 2. Integrating up to z = 10 one 
can determine AG(fi 2 ) with 20% accuracy for the distribution of Brodsky et al.|l2|, while 
for the distribution of Chiappetta et al. |L3| the accuracy is even better, of the order 



of 10%. Note that the distribution of Ref. [13| yields a Ioffe-time distribution T(z) which 
oscillates around zero at large z. Therefore a fraction of the full integral (|9]) obtained by 
an integration up to a finite value of z can be larger then one. The slightly worse accuracy 
obtained for the distribution of Ref. W2 is due to the fact that in this case the polarization 
AG([i 2 ) is much smaller and the corresponding Ioffe-time distribution has a rather flat 
shape, which is seen when both distributions are compared at the same vertical scale. 

We show now that only a few moments of AG(x, /j, 2 ) are necessary to obtain a rea- 
sonable estimate of AG(/i 2 ). Already the first two non- vanishing moments are sufficient 

1 Note that the QCD interpretation of is similar to the Gottfried sum rule which also cannot be 
connected to a matrix element of a twist-2 local operator, but is determined by an infinite number of 
matrix elements of a tower of twist-2 operators. 

2 Such a picture makes sense only as far as the nucleon Compton wavelength is about 5 times shorter 
than its electromagnetic radius. 
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to obtain the approximation to the Ioffe-time distribution shown in Figure 1 and 2: 

1 



6 



(10) 



As an approximation to AG(/i 2 ) we calculate the area bound by the approximate Ioffe- 
time distribution (10) up to its maximum, and the straight line connecting this point with 
the value of T(z = 10, /i 2 ). For the distributions of Ref. [112] and [Tj| this procedure yields 
AG = 0.41 and AG = 1.67, as compared to the exact values 0.54 and 1.7, respectively. 
Although due to the reasons discussed above the accuracy is lower in the first case, in 
general both estimates are quite satisfactory. In an even simpler approach, calculating 
the area of the triangle spanned by the points z — 0, z — 10, and the maximum of the 
approximate Ioffe-time distribution (|i~0"|), one obtains: 



ag( m 2 ) » ^r 2 (^) 



2r 2 (/i 2 



'ID 



which yields AG = 0.28 for the distribution of Ref.jl^] while for Ref.[[L3 
AG = 1.1, i.e. an accuracy of about 50% in both cases. 



one obtains 
This demonstrates that the 
gluon polarization is indeed dominated by the two first moments of the polarized gluon 
distribution. We note that the whole procedure assumes implicitly that the Ioffe-time 
polarized gluon distribution has a regular shape which in principle cannot be taken for 
granted. On the other hand all presently advocated distributions |14| have this property, 
and therefore it can be considered at least as a reasonable hypothesis. 

As explained above, besides (at least) the first two moments, we need to know the 
normalization of T(z, /i 2 ) at large z ~ 10. As the large z region corresponds to the small x 
domain, we need information about the shape and normalization of the small x polarized 
gluon distribution. There, the color coherence of gluon couplings results in a simple 
relation between the unpolarized and polarized gluon distribution 



'AG(x) 



■r. 



for x — > 0. 



nucleon 



(12) 



Using GRV |15| and CTEQ |L6| LO unpolarized gluon distributions we have found that 
T(z = 10) = 0.005 - 0.007 at 1 GeV 2 . 

As far as moments of twist-2 parton distributions are concerned, lattice calculations 
have produced already quite satisfactory predictions for the first two moments of unpo- 
larized and polarized G 



+ 1 and G = — 1 quark distributions IT? 



matrix elements of gluonic operators are not yet as reliable |18 



but extensions to 
Hence, to illustrate the 



above procedure we have computed the second and fourth moment of the polarized gluon 
distribution at a scale /i 2 ~ 1 GeV 2 using QCD sum rules 
calculation of matrix elements of QCD operators |2T) |. 



19], suitably tailored for the 



For this purpose we consider the three-point correlation function 



'AG 



d A xe^ x / d*ye i ™{0\T[r ]G (x)O(y + A/2;y - A/2)t? g (0)]|0) (13) 
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of the operator (|5]) with isospin 1/2 interpolating currents T]g(x) which contain explicitly 
gluonic degrees of freedom: 



r] G {x) 



^(Vo\d(x) -7] ex (x)), 



^abc/„.aT 



r abc( „,aT 



u aT {x)C 1 ^d b {x)) lh ^a aP \gG^{x)u{x 



(14) 



The current t]g{x) was first introduced in |^TJ to compute the first moment of the unpo- 
larized gluon density. Since then it has been applied to estimate higher twist corrections 



to the Bjorken sum rule |p2| , and to calculate the nucleon matrix element of the pseudo- 
scalar gluon density ^-G 11V G^ U |23|, giving in both cases reasonable results. From a 



technical point of view the main advantage in using the current t]g is that one avoids 
UV divergences which occur when gluons are generated perturbatively and which survive 
the Borel transformation. Although in principle we are interested only in the first two 
local operators which contribute to the Taylor expansion of (|5|) around A = 0, it turns 
out that one can account for the non-local operator 0(y + A/2; y — A/2) as a whole at 
no additional expense. To avoid significant non-perturbative contributions from large t- 
channel distances we have simply kept Q 2 = —q 2 ~ 1 — 3 GeV 2 in the Euclidean domain, 
and extrapolate numerically at the end to Q 2 = 0. The kinematics is chosen such that 
q • A = and q 2 = —q±, which yields p ■ A = {p + q) ■ A = z. 

We focus on the contribution to the correlator ([TBI) of the form: 



(15) 



lbP{p-n) 2 [ duT(p 2 ,(p + q) 2 ,Q 2 ,u)sin(uz) 
Jo 



where p = p^Tf^ 1 . The invariant function T(p 2 , (p+q) 2 , Q 2 , u) can be projected out uniquely 
by taking the trace of Jag with jfij^. Identifying the contribution from nucleon interme- 
diate states to the three-point correlation function (^) leads to: 



l Tr (w7 5 /ag) = I duT(p 2 ,(p + q) 2 ,Q 2 ,u)sm(uz), 
4 (p ■ n) A Jo 

-T(z,Q 2 ) + 



\ 2 G m% 



N 



p 2 )(m 



n - (p + q¥ 



(16) 



where the form-factor T(z;Q 2 ) coincides at Q 2 = with the Ioffe-time distribution T(z) 
in eq.(|^). In ( Jl6| ) the nucleon polarization has been fixed parallel to the 3-axis such that 
s ■ A = p ■ A, and ellipses denote contributions from higher resonances and the continuum. 
Note that here and in the following we have suppressed the explicit reference to the 
renormalization scale fi 2 which is understood to be approximately equal to 1 GeV 2 . The 
overlap of the interpolating current i]g with the nucleon state, 



(0\r) G \p,s) 



m 



N 



X G u{p, 



(17) 
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at the scale /i 2 ~ 1 GeV 2 has been determined in Ref. pi 



The invariant function T(p 2 , (p + q) 2 ,Q 2 ,u) admits a double spectral representation 



2 r dsi f ds 2 



T(p\ (p + q)\ Q, u) = / L — / -^—^ p( Sl , s 2 , Q\ u). (18) 

Si - p 2 J s 2 - {p + q) 2 



Therefore, after expanding the left- and right-hand side of eq . ([lq) around A = we 
obtain: 

T l{ p\ {P + q ?^) = J^J-^—^^ 

= r~2 2^7 % ( 7 , 2 M Q 2 ) + ..., / = 2,4.... (19) 

(m 2 N - p 2 ) (m 2 N - (p + <?) 2 ) 

Here p/(si, S2, Q 2 ) = J duu 1 ^ 1 p(si, s 2 , Q 2 , u) is the /-th moment of the spectral density. A 
comparison with eq.(|4]) shows that in the limit Q 2 — > the moments Ti(Q 2 ) are identical 
to the moments of the polarized gluon distribution. To eliminate the contributions from 
higher resonances and the continuum we limit, according to the standard procedure, the 
integrals over s± and s 2 by the continuum threshold sq and subsequently perform a Borel 
transformation in both p 2 and (p + q) 2 , arriving at the sum rule: 



\ 2 G m% 

Note that the continuum threshold sq and the Borel parameter M 2 should in principle be 



UQ 2 ) = "T^T" / ds ± / ds 2 e" (si+S2)/2M2 pK^i, s 2 , Q 2 ) , (20) 
\ r m\j Jo Jo 



kept at their values fixed by the two-point nucleon sum rule [21 



In the following we concentrate on the sum rules for the first two moments of AG(u). 
In the calculation of the spectral function p(si, s 2 , Q 2 , u) we have restricted ourselves, as 
in plf, to the perturbative and the four-quark condensate contributions only. Explicitly, 



we obtain: 

rf"V.*.«") = 28 °^ fl^ (5fa + ^)fl-22 SlS ^), (21) 
= ^( 5 ,)^£l£^!( fl + 7. SlS2 ), (22) 

pf\si,S2,Q 2 ) = ^ (i'j) 2 "^tf' 2 (fi 2 + 18^i«2 + <»44) ■ (23) 

where A q = si + s 2 +Q 2 and R = A 2 — 4siS 2 . The expression for P2 ert ( s i> s 2, Q 2 ) is unfortu- 
nately too long to be quoted here. It is interesting to note that in this 
approximation the calculation of the three-point function with the operator 
7i M n^TrG^(A/2)[A/2;-A/2]G£(-A/2), determining the Ioffe-time distribution of un- 
polarized gluons, results in exactly the same spectral density p(si, s 2 , Q 2 , u). In numerical 
calculations we have used —(27i) 2 (qq) ~ 0.67 GeV 3 which is the standard value of the 
quark condensate at a scale of 1 GeV 2 . 
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1/2 rp i/i 

s [GeVJ 


1V1 [Ue V J 


J- 2 


J- 4 


1.4 


1 


0.289 


0.205 


1.5 


1 


0.340 


0.240 


1.6 


1 


0.387 


0.265 


1.4 


2 


0.312 


0.217 


1.5 


2 


0.389 


0.273 


1.6 


2 


0.477 


0.335 



Table 1: T 2 and T 4 for various values of sum rule parameters. 



To determine Y\ = Ti(Q 2 = 0) we have fitted I^Q 2 ) and T^(Q 2 ) in the interval 
1 < Q 2 < 3 GeV 2 with the ansatz: 

r > (Q2) = WTWW>' (24) 

as suggested by spectator counting rules. To find the dependence of our results on s and 
M 2 we have varied them around the standard values known from the two-point nucleon 
sum rule. The estimates for r 2 and r 4 obtained in this way are collected in Table 1. We 
observe a moderate dependence on the sum rule parameters. Taking ^/sq = 1.5 GeV and 
M 2 = 1 GeV 2 we get T 2 = 0.34 and T 4 = 0.24 at the scale /i 2 ~ 1 GeV 2 . 

With the above results at hand, we obtain the Ioffe-time distribution T(z, fx 2 ~ lGeV 2 ) 
in the small z domain by the expansion in eq. fliCf) . Together with our estimate for T{z — 
10) discussed above, we get AG « 2.0 ±0.5 at the scale /i 2 ~ 1 GeV 2 . The error estimate 
results from the quality of the sum rule. This value is not far from AG(lGeV 2 ) = 1.7 
for the parameterization of Chiappetta et al. |13[ , and it supports the recent analysis p4 



of the Q 2 -dependence of polarized structure function data which yields AG(1 GeV 2 ) = 
1.5 ± 0.8. Contributions to T(z) from z > 10 can change the value of AG by about 
10%. On the other hand the present calculation overestimates the true values of T/. This 
can be seen for example from our prediction for T 4 being larger then the experimentally 



determined value of the third moment of the unpolarized gluon distribution f25| , which 
contradicts unitarity. Our result for AG should therefore be interpreted as an upper limit. 

To summarize, we have found that for sufficiently smooth polarized gluon distributions 
the gluon polarization in the nucleon can be estimated to good accuracy from the first 
two even moments of the polarized gluon density, together with its behavior at large Ioffe- 
times, or equivalently at small values of Bjorken variable x. We have calculated the former 
within the framework of QCD sum rules. Taking into account the perturbative and the 
four-quark condensate contribution we found the corresponding spectral function to be 
similar to the spectral function for the gluonic operator determining the unpolarized gluon 
distribution. Consequently, since the contribution of gluons to the nucleon momentum 
is around 50% even at low scales /i 2 ~ 1 GeV 2 , also the polarized gluon distribution is 
expected to be large. Combining the calculated moments with the color coherence of 
gluon couplings we have obtained AG(fi 2 ~ 1 GeV 2 ) ~ 2, which is acceptable from a 
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phenomenological point of view. Finally we want to emphasize the need for a calculation 
of the first two moments of the polarized gluon density on the lattice. Unfortunately, 
even an order of magnitude estimate of T 4 requires the treatment of a symmetric traceless 
tensor of rank 5 and therefore appears to be quite difficult. 
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Fig. 1 The solid line represents the Ioffe-time distribution T(z) of the polarized gluon 
density of Brodsky et al. The approximate Ioffe-time distribution defined by 

the moments and T4 flit] ) is shown by the dotted line. The fraction of the first 
moment integral @ as a function of the upper integration limit z is given by the 
dashed-dotted line. The shaded area yields the approximation to AG as explained 
in the text. 
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